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We reanalyze the cosmological constraints on the ex- 
istence of a net universal lepton asymmetry and neu- 
trino degeneracy based upon the latest high resolution 
CMB sky maps from BOOMERANG, DASI, and MAXIMA- 
1. We generate likelihood functions by marginalizing over 
(fl^h 2 , £„ T , £„ e , Q,a, h, n) plus the calibaration uncertainties. 
We consider flat Qm + = 1 cosmological models with two 
identical degenerate neutrino species, £i/ M>T = | = \£,v T 
and a small £„ e . We assign weak top-hat priors on the 
electron-neutrino degeneracy parameter £„ e and f2j,/i 2 based 
upon allowed values consistent with the nucleosynthesis con- 
straints as a function of £„ MjT . • The change in the background 
neutrino temperature with degeneracy is also explicitly in- 
cluded, and Gaussian priors for h = 0.72±0.08 and the exper- 
imental calibration uncertainties are adopted. The marginal- 
ized likelihood functions show a slight (0.5<r) preference for 
neutrino degeneracy. Optimum values with two equally de- 
generate fi and t neutrinos imply = 1.0^°'*^ , from 

which we deduce £„ c = 0mt°m, and Q b h 2 = 0.021±g'go2- 
The 2a upper limit becomes T < 2.1, which implies 
£ Ve < 0.30, and Qth 2 < 0.030. For only a single large- 
degeneracy species the optimal value is \£v \ or |£„ T | = 1.4 
with a 2a upper limit of |£„ | or |£„ T | < 2.5 



I. INTRODUCTION 

The present relic neutrino number asymmetry is not 
directly observable. Hence, there is no firm experimental 
basis for postulating that the lepton number for each 
species is zero. Charge neutrality, however, demands 
that any universal net lepton number beyond the net 
baryon number must reside entirely within the neutrino 
sector. It has been suggested that the total lepton num- 
ber could be large in the context of the SU (5) and 50(10) 
grand unified theories , or supersymmetric baryoge- 
nesis ^j-Q based upon the Affleck-Dine scenario || . Such 
mechanisms might generate lepton number asymmetry 
up to ten orders of magnitude larger than the baryon 
number asymmetry. Furthermore, even if one demands 
that B — L « 0, it is possible for the lepton numbers 
Li of individual neutrino species to be large compared to 
the baryon number of the universe, B, as long as the net 
total lepton number is small. 



Moreover, there presently exists at least some marginal 
cosmological evidence for neutrino asymmetry. For 
example, neutrinos with large lepton asymmetry and 
masses ~ 0.07 eV might be required to explain the ex- 
istence of cosmic rays with energies in excess of the 
Greisen-Zatsepin-Kuzmin cutoff [§ 10|- Also, degenerate, 
massive (2.4 eV) neutrinos might be required |11| to pro- 
vide a good fit to the power spectrum of large scale struc- 
ture in mixed dark matter models. It is thus important 
to carefully scrutinize the limits which cosmology places 
on the allowed values of possible universal lepton asym- 
metry. Indeed, a number of recent papers [fl2|-p3f have 
addressed this issue with varying degrees of complexity. 
The present work differs from those in several details as 
summarized below. It represents an independent exami- 
nation of this issue. 



II. PRESENT APPROACH 

In a recent paper |24j ] we considered new constraints 
imposed on neutrino degeneracy from primordial nucle- 
osynthesis. Particular attention was paid to the neu- 
trino decoupling temperatures before the nucleosynthesis 
epoch. Of relevance to the present work is that we have 
shown that neutrinos can decouple at a higher temper- 
ature than estimated in earlier studies ]25| ]. This means 
that more particle degrees of freedom could be present at 
neutrino decoupling. This causes the relic neutrino tem- 
perature to be lower by simple entropy considerations. A 
smaller relic neutrino energy density implies that larger 
neutrino degeneracies may be allowed. For example, we 
have shown that interesting regions of the model parame- 
ters for big-bang nucleosynthesis (BBN) are allowed such 
that substantial lepton asymmetry and baryon density 
(even ttbh 2 w 1 where h is the present value of the Hub- 
ble constant in units of 100 km sec -1 Mpc -1 ) are possible 
while still satisfying the adopted abundance constraints 
from primordial nucleosynthesis. 

A stronger constraint, however, on lepton asymmetry 
comes from the power spectrum of fluctuations in the 
cosmic microwave background (CMB) which we now ad- 
dress in the present paper. We apply a likelihood analy- 
sis of neutrino-degenerate models to the combined latest 
BOOMERANG [§|, DASI Q and MAXIMA-1 |2§] re- 
sults. We note, however, that a recent analysis p9] of 
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the implications of neutrino oscillations derived from a 
combination of the atmospheric, and solar neutrino con- 
straints implies much tighter limits on degeneracy for all 
neutrino flavors. If neutrino oscillation parameters are in 
the range of the large mixing angle solution then an up- 
per limit of 1^0.07 applies to all neutrino flavors. The 
limits derived here do not assume any particular model 
for neutrino mixing, and should be taken as independent 
of, and complementary to, those constraints. 

The implications of the CMB data for neutrino- 
degenerate cosmologies have been noted in a number of 
recent papers ]l6|-p3|]. The constraints on the effective 
number of relativistic particles can also arise in other con- 
texts, such as cosmic quintessence |30],[3l|. The present 
work, however, differs from those in several respects. For 
one, we consider the most recent combined data sets, 
not just the first year BOOMERANG data as in ff|[L7). 
and generate a marginalized likelihood function for the 
neutrino degeneracy and other cosmological parameters. 
Many of the existing studies have marginalized over a 
more limited set of cosmological parameters. For exam- 
ple, in ]l9| only fi^h 2 and and neutrino degeneracy were 
marginalized to set limits while other cosmological pa- 
rameters were set to various fixed values. In [V7\, for 
example, no likelihood analysis was made. In [jl6| a like- 
lihood analysis was made but without window functions. 

For the present work we use the updated Radpack 
package |3^] described below which includes all relevant 
window functions. Another difference between our anal- 
ysis and other works is that our marginalization utilizes a 
global minimum search algorithm ]33[ ] rather than a dis- 
crete grid of cosmological parameters. Marginalizing pa- 
rameters for each fixed value of one parameter requires at 
least 1000 model calculations to get 10 -4 accuracy for the 
X 2 minima even by using this algorithm. Nevertheless, 
in this way we are sure to identify the true marginalized 
likelihood functions. 

The most similar recent likelihood analysis to that de- 
scribed here is in the work of [|T). Our analysis differs 
from [pl|] in several respects. In the present work we 
make use of our deduced new family of baryon densi- 
ties and lepton asymmetries allowed by BBN to assign 
weak top-hat priors on the derived likelihood functions. 
This differs from that in |pl[| in which separate Gaussian 
likelihood functions were evaluated for the nucleosynthe- 
sis constraints and the CMB. A total likelihood function 
was then defined by marginalizing over the product of 
these two functions. We prefer our method because the 
uncertainties in the BBN constraints are dominated by 
systematic errors. Systematic errors are not equivalent 
to random Gaussian errors. We thus prefer weak top- hat 
priors as a more realistic representation of the systematic 
errors in the BBN constraints. 

One other important difference is that we adopt a 
strong Gaussian prior of h = 0.72 ± 0.08 based upon the 
Hubble Key Project results Q. In [plj weak a top-hat 



prior of h — 0.65 ± 0.20 was adopted. As noted above, 
another difference between the present work and all pre- 
vious results is that we consider carefully the change in 
background neutrino temperatures as a function of de- 
generacy. Although this is a small effect for low degen- 
eracies, it can slightly affect the upper limits. 



III. NEUTRINO-DEGENERATE BBN 

Neutrino degeneracy affects BBN in two ways. The in- 
clusion of a small amount of electron v e degeneracy resets 
the equilibrium neutron to proton ratio at weak-reaction 
freezeout to n/p — exp{— Am/T — £,v,,}- This can cause 
a reduction in the primordial helium abundance. Indeed, 
it has been argued J|5|,[56| that the apparent conflict be- 
tween the low helium abundance inferred from HII re- 
gions of metal poor galaxies and the low Lyman-a deu- 
terium abundance may even require v & degeneracy for its 
resolution. The present deuterium-absorption limits on 
ilb/i 2 « 0.020 ±0.002 (2a) requires a large primordial he- 
lium abundance of Y p ^0.25 and substantial destruction 
of primordial 7 Li in stars. These conditions tax even 
the most generous adopted limits from observed light- 
element abundances 37|. Thus, a modification of BBN 
which allows for large values of Vt^h 2 while still satis- 
fying the constraints from light-clement abundances is 
worth investigating. Such conditions are easily satisfied 
by neutrino-degenerate models. 

The inclusion of either or v T degeneracy on the other 
hand only enhances the background energy density and 
therefore the universal expansion rate. During the radi- 
ation dominated epoch, relativistic neutrinos contribute 
a large fraction of the mass energy. Thus, even a small 
modification of the neutrino energy density can signifi- 
cantly affect the expansion. 

The energy density p v due to degenerate neutrinos (or 
any other fermions) are described by the usual Fermi- 
Dirac distribution functions /„ = [exp (EjTy — £„) + 
where the neutrino degeneracy parameter is defined 
by ^ v = ii v /T v , and [i v is the neutrino chemical potential. 
Thus, we have 



Pp = 
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d PP 2 E u (fAp) + U(p))- 



(1) 



where, p denotes the magnitude of the 3-momentum, and 
E v = \J p 2 + m 2 , with m v the neutrino mass. Here and 
throughout the paper we use natural units ( h = c = ks = 
I)- 

For the present discussion it is sufficient to only con- 
sider massless neutrinos. (Possible limits on neutrino- 
degenerate models with massive neutrinos are considered 
in [p^[). The energy density in massless neutrinos be- 
comes 
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from which it is clear that degeneracy in any neutrino 
species tends to increase the energy density. The associ- 
ated increased expansion rate tends to increase the neu- 
trino decoupling temperature. This causes an increase 
in the primordial helium and other light-clement abun- 
dances. 

Since £„ and f„ T primarily affect the expansion rate, 
they are roughly interchangeable as far as their effects on 
nucleosynthesis or the CMB are concerned. Furthermore, 
it now seems likely [p8[ that the mixing parameters for 
and v T involve a large mixing angle and small 8m 2 . 
In this case it is plausible that the muon and tau neutri- 
nos were interconverted in the early universe and would 
therefore obtain nearly an identical degeneracy parame- 
ter if an asymmetry exists. Thus, we adopt a conserva- 
tive model in which the \i and r neutrinos are equally 
degenerate, %^ | = |£„ T | = . 

As shown in |24j , for each value of T there is a 
unique range of £„ e and Sl^h 2 which satisfies the com- 
bined deuterium and primordial helium constraints. The 
allowed family of neutrino-degenerate models employed 
in this work is summarized in Figure |Tj. 

This figure differs slightly from the family of allowed 
solutions given in [ p4| in that we have adopted the newer 
D/H constraint from |$| (i.e. D/H = 3.0(±0.4) x lfr 5 ) 
and slightly different limits on the primordial helium 
abundance (0.228 < Y p < 0.248). In the limit of the 
standard nondegenerate big bang (£„ = £j, e = 0) our 
limits on fl^ti 2 reduce to those of [{40| , f4~L|| , i.e. fi^h 2 = 
0.021 ± 0.002. The allowed shaded regions in Figure |l| 
will be adopted as weak top-hat priors in the CMB like- 
lihood analysis described below. These regions include 
both the uncertainties from the abundance constraints 
described above and the uncertainties in the BBN model 
predictions p2|. 



IV. CMB POWER SPECTRUM 

Having defined the family of allowed priors from BBN 
we can now do a likelihood search for optimum cosmolog- 
ical parameters which fit the CMB data. Several recent 
works ]l2]-|l6| have explained how neutrino degeneracy 
can dramatically alter the power spectrum of the CMB. 
For massless neutrinos it can be shown |l5|] that the only 
effect of neutrino degeneracy is to increase the back- 
ground pressure and energy density of relativistic par- 
ticles. The essence of this constraint is that degenerate 
neutrinos increase the energy density in radiation at the 
time of photon decoupling and delay the time of matter- 
radiation energy-density equality. This mainly causes an 



increase in the amplitude of the first acoustic peak in the 
CMB power spectrum at I « 200. For example, based 
upon a \ 2 analysis Jl3| of 19 experimental points and 
window functions, it was concluded in Jill that £„ < 6 
for a single degenerate neutrino species with an Q\ = 
cosmology. 

However, in the existing CMB constraint calculations 
]l2| p3| only small degeneracy parameters with the stan- 
dard relic neutrino temperatures were studied in the 
derived constraint. Hence, the possible effect of a di- 
minished relic neutrino temperature at high degeneracy 
needs to be considered. To investigate this we have 
done calculations of the CMB power spectrum, AT 2 = 
l(l+l)Ci/2ir based upon the CMBFAST code of Seljak & 
Zaldarriago . We have explicitly modified this code to 
account for the contribution of massless degenerate neu- 
trinos with varying relic neutrino temperatures T Vi for 
each species Q. 

The experimental uncertainties are non-Gaussian, but 
can be well represented by an offset log-normal distribu- 
tion [||| . As in we have evaluated the x 2 goodness of 
fit for a range of theoretical power spectra C; as follows: 
We define the goodness of fit by 



^(Zt-Zf)Mjj{Zl-Zf) 



Xcal 



(3) 



where separate summation over the different data sets is 
implied. For each set of binned power data Zf we utilize 
the published off-set log-normal data from the three data 
sets. 
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Where Di is the measured band power. The correspond- 
ing binned theoretical power spectra are 



Z* = ln 



Y,^(Wu/l)(C l + x l ) 



(5) 



where the are the published calibration uncertainties 
taken to be 8%, for MAXIMA-1 and DASI, and 20% 
for BOOMERANG. Window functions W u for the three 
data sets are available on the world wide web. The error 
matrix is simply 



Ml 



MijiPi+x^iDj+Xj) 



(0) 



where Mij is the weight matrix for the band powers D^. 
The effect of the calibration uncertainty on the goodness 
of fit is obtained from 
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where <7j is the experimental uncertainty. The total \ 2 
evaluated in this way can be converted (27].^(| into a like- 
lihood function for each parameter x marginalized over 
the remaining parameter set y 
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C(x)= P P rior(x,y)exp(-x 2 /2)dy 



(8) 



In neutrino-degenerate models which satisfy the con- 
straints from primordial nucleosynthesis |2~i| ], increasing 
the neutrino-degeneracy must be accompanied by a com- 
mensurate increase in baryon density. Fits to the CMB 
power spectrum for large degeneracy |24| , therefore show 
a suppression of the the second acoustic peak due to 
baryon drag [jl2| . 

Indeed, such suppression of the second acoustic peak 
seemed to be present in the first reported power spec- 
tra based upon the balloon-based CMB sky maps from 
the BOOMERANG @ and MAXIMA-1 @] collab- 
orations. This remains true for the likelihood anal- 
ysis based upon a MAXIMA-1 data |^5| which indi- 
cates fl b h 2 = OMOtoQw (2cr). However, in the most 
recent data sets from' BOOMERANG || and DASI 
f27f the second peak has become much better defined. 
Both the BOOMERANG and DASI data sets now im- 
ply n b h 2 = 0.022l°;<$3 (Iff) (7710 = 6.00l£i°). This 
value is close to the value implied by the cosmic deu- 
terium abundance in high-redshift Lyman-a clouds ob- 
served along the line of sight to background quasars 
@,|3 n b h 2 = 0.020 ± 0.001 (lcr) (7710 = 5.46 ± 0.27). 
Hence, the newer data imply at most a marginal require- 
ment for a larger baryon density or neutrino degener- 
acy. Indeed, these new data tighten constraints on the 
possibility of degenerate cosmological neutrinos. In the 
present paper we explore the new limits on possible neu- 
trino degeneracy implied by the combined data sets and 
our BBN constraints. 



V. RESULTS 

We limit our consideration to flat £ltot — &m + = 1 
cosmological models with ionization parameter r = 0. 
This is sufficient for our purposes since the likelihood 
functions so deduced are not expected to be much dif- 
ferent if Vitot or r are varied (cf. 26 1). This is because 
r and the spectrum tilt n are nearly degenerate parame- 
ters, i.e. changing one is equivalent to changing the other. 
Moreover, f2 tot is generally tightly constrained to be near 
unity anyway. 

There are then nine parameters over which we 
marginalize. These are (fl b h 2 , £„ M T , ^ c ^a, h, n, e^). We 
utilize a strong Gaussian prior for h = 0.72 ± 0.08 and 
for the calibration uncertainties Ci as listed above. Also, 
as noted above, we adopt weak top-hat priors when 
marginalizing over fl b h 2 and £„ e designated by the shaded 
regions of Figure |l| for each value of £ v T . In |2^] it has 
been argued that without some priors on Qm (through 
flatness, h, etc.) it is difficult to place bounds on the 
amount of relativistic matter. Hence, the model con- 
straints adopted here are probably required to break the 



parameter degeneracy between relativistic and nonrela- 
tivistic matter. Ultimately, however, high resolution sky 
maps such as the Planck mission will be able to determine 
separately the amounts of relativistic and nonrelativistic 
matter. 

Figure [2] illustrates one of the main results of this 
study. Shown are contours of constant A% 2 in the £ y 
vs. n plane for three values of the cosmological constant 
(i.e. n A = 0.65, 0.75, and 0.8) and for fixed h = 0.75 
as noted. For J7a < 0.75, a minimum in x 2 develops for 
values of T w 1 — 1.5. Indeed, for a simple 2 parameter 
search with fixed values, f^A = 0.75 and h — 0.75, neu- 
trino degeneracy is preferred at the level of more than 3<r 
over a nondegenerate model. For smaller values of f^A, 
this minimum for neutrino-degenerate models becomes 
even more pronounced. 

A second minimum also develops for higher degeneracy 
(6^ ~ 11.4) as noted in [pi} . This is due to the large 
change in particle degrees of freedom for neutrinos which 
decouple just above the QCD transition. However, the 
goodness of fit is so poor (A% 2 > 500) that it would not 
be apparent in the contours drawn on Figure |[ Hence, 
this large-degeneracy solution is definitely ruled out by 
the current CMB power spectrum. 

Figure |^ shows the marginalized likelihood dis- 
tributions for three of the cosmological parameters 
{£ v , f^A, n) considered here. For the present study 
the likelihood functions for £„ e , Q, b h 2 and VLm are related 
to these since £ Ve and fi b h 2 are functions of £ v T and 
Qm = 1 — ^a- From Figure || we deduce optimum values 
of n A = 0.74l£?i and n = 0.93 ± 0.02. A slight pref- 
erence for finite neutrino degeneracy is evident T = 

l-O^i o(o5 ff )' This preference, however, is not particu- 
larly significant. For now, the data mainly imply (2a) 
upper limits on neutrino degeneracy of £„ < 2.1. This 
value implies upper limits of £„ e < 0.30 and fl^h 2 < 0.030 
from Figure [l]. For a single large-degeneracy neutrino 
species, these limits become |£j,J or |£„J = 4/J 
with a 2a upper limit of £„ < 2.5. 

Our results are slightly more stringent than the results 
from |2~i|| who found an equivalent single species upper 
limit based upon the CMB data alone of £„ or £ Vt < 2.9. 
This is at first surprising given that we have adopted 
weak top-hat (instead of Gaussian) priors for the BBN 
constraint. We have traced the main reason for the more 
stringent upper limits derived here to our adoption of a 
strong Gaussian prior on h. A larger neutrino degeneracy 
is possible if larger values of h are permitted. Figure |] 
shows contours of constant A% 2 in the Ho vs. n plane for 
f?A = 0.75 models with = 0, 1.0 and 1.5 as labeled. 
This illustrates the sensitivity of the degenerate solution 
to the assumed prior for h. If weaker priors on h are 
adopted, or if new larger values of h in the upper range of 
the present Key-Project uncertainty are ever determined, 
the neutrino-degenerate models could become strongly 
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preferred over the non-degenerate models. The = 
non-degenerate solution is only the preferred minimum, 
for all values of fi/v, when h < 0.70. This is consistent 
with the results of [||,|6). 

Figure [s] shows some optimum model power spectra 
compared with the combined data sets. The solid line 
shows our optimum degenerate model for which (fl^h 2 , 
&,„, T , n A , h, n) = (0.021, 1.0, 0.09, 0.74, 0.74, 0.93). 
For this parameter set we obtain a total \ 2 — 29.8 for 29 
degrees of freedom implying a nearly perfect fit. For com- 
parison, the dotted line shows the best non-degenerate 
(& MlT = Zv t = 0) model [{VL b h 2 , n Aj h, n) = (0.021, 0.62, 
0.62,1.0) (dotted line)] from pq| . For illustration we also 
show the large-degeneracy minimum [(il^/i 2 , £„ , £„ e , 
n Aj /i, n ) = (0.052, 11.4, 0.74, 0.45, 0.80, 0.72) (dot- 
dashed line)]. 



VI. CONCLUSIONS 

In neutrino-degenerate models the larger baryon den- 
sity associated with the observed low deuterium abun- 
dance can be more easily accommodated than in non- 
degenerate models. Moreover, neutrino-degenerate mod- 
els provide a slightly improved goodness of fit for the 
latest CMB power spectra from BOOMERANG, DASI, 
and MAXIMA- 1. 

Using cosmological models consistent with the con- 
straints from light-element abundances as a function of 
the neutrino degeneracy parameter T , we have shown 
that a slight maximum in the likelihood function forms 
for neutrino-degenerate models with w 1. However, 
the improvement over the nondegenerate models is only 
at the level of about 0.5er. Although this minimum is not 
particularly statistically significant for the present data 
set and assumed priors, it could become much more pro- 
nounced should larger values of h and/or smaller values 
of Q\ ever be established near their current la limits. 

The present data place 2a limits for two identical large- 
degeneracy neutrino species is £„ < 2.1, which implies 
£i/ e < 0.30. For only one species with large degeneracy, 
the limit becomes \£ v | or 1^ | < 2.5. This is slightly 
more restrictive than the limits deduced in other studies. 

Finally, we remark that, since neutrino degeneracy is 
now limited to such small values, the present work has 
established that the effects of the changing neutrino de- 
coupling temperature with increased degeneracy has lit- 
tle effect. Hence, previous studies which neglected this 
effect are justified. 
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FIG. 1. Allowed values of Q b h 2 and f„ e for which the con- 
straints from light-element abundances are satisfied as a func- 
tion of £„ . 
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FIG. 4. Contours of constant goodness of fit A\ 2 in the 
Hq vs. n plane for Qa = 0.75 and a neutrino degeneracy 
parameters £i/ MtT . = 0., 1.0, and 1.5 as labeled. 
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FIG. 5. Fits to the power spectrum of fluctuations in the 
CMB. The solid line shows the best neutrino-degenerate fit 
t — 1-0). The dotted line shows a best non-degenerate 
t — £v e = 0) model. For illustration, the dot-dashed line 
also shows the large-degeneracy minimum (f„ T = 11.4). 



8 



